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I. INTRODUCTION 

Since the pioneering work of Coleman and De Luccia (CdL) [1] many articles were devoted 
to the investigation of metastable vacuum decay with gravity. The recent discovery of a Higgs 
boson at the LHC highlights the importance of knowing all possible corrections (including 
gravitational ones) to the vacuum decay rate, given that the present Higgs boson and top 
quark mass values indicate that we may live in a metastable vacuum [2-14], In spite of many 
investigations of various aspects of metastable vacuum decay with gravity^, there remain 
several important open questions in the held. One of them is the negative mode problem. 
While in hat space-time it is relatively easy to show that a bounce solution describing 
metastable vacuum decay [19-21] contains exactly one eigenmode with negative eigenvalue 
(we refer to this as the tunnelling negative mode) in its spectrum of linear perturbations 
(as it should - see [22]), once gravity is included this question is still not completely solved 
[23-37]. 

The quadratic action for small perturbations about CdL bounces was hrst calculated 
by Lavrelashvili, Rubakov and Tinyakov (LRT) [23] using a Lagrangian formulation and 
hxing the gauge in order to obtain an action for the single remaining degree of freedom. It 
was found that this action has the following structure: one obtains a factor Qlrt, given 
by Eq. (26) below, in front of the kinetic term of the perturbations. This factor Qlrt 
depends on a certain combination of background quantities and typically becomes negative 
on some interval along the bounce trajectories. This specihc gauge choice was criticised in 
[24], where it was argued that the gauge is ill-dehned at the middle of the instanton, i.e. at 
the location of maximal 3-volume of the geometry.^ In fact, in the Lagrangian approach, 
where only conhguration variables are involved, many gauge choices lead to similar kinds of 
trouble. Motivated by this, Tanaka and Sasaki (TS) [24] derived a quadratic action using the 
Hamiltonian formalism. After excluding matter degrees of freedom TS obtained the action 
for one (gravitational) degree of freedom and showed that there are no negative modes 
associated with this quadratic action [24, 26]. Realizing that the approach of TS is not fully 
satisfactory, in particular because it does not allow one to recover the flat space limit when 

^ Aspects of interest, but not dealt with in detail here, include [15-18]. 

^ This is reminiscent of the issue of gauge fixing in bouncing universes, where the scale factor also passes 
through an extremum - see e.g. [38]. 
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gravity is turned off and because it hinges on a number of delicate analytic continuations, a 
new approach was suggested by Khvedelidze, Lavrelashvili and Tanaka (KLT) [27]. Using the 
full machinery of constrained dynamical systems, the quadratic action for the single physical 
degree of freedom was obtained, and in this framework the flat space limit is successfully 
recovered. Note that similar results were later re-derived in [29] and [32]. The end result 
of the KLT approach is similar to the earlier Lagrangian approach action, but it contains 
a different factor Q in front of the kinetic term of the fluctuations, see Eq. (27) below. It 
has been proven [27] that for bounces where Q is everywhere positive there exists exactly 
one negative mode in the spectrum of linear perturbations (see also [28, 29, 32]). In this 
case, in complete analogy with the non-gravitational case, one can interpret the instanton 
as mediating the decay of a metastable vacuum. 

A recent study of gravitational instantons in flat potential barriers has highlighted the 
fact that for large classes of instantons, Q becomes negative somewhere along the instanton 
[35]. This has motivated us to study the issue of negative modes in some detail for these 
cases. Our findings are twofold: 

First, we will demonstrate (analytically) that the eigenvalue equation retains non-singular 
solutions when Q passes through zero. In other words, at Q = 0 the equation is seemingly 
singular, but its solutions remain continuous and at least twice differentiable. Moreover, the 
eigenvalue evolves continuously as the potential is varied from cases where Q > 0 everywhere 
to cases where Q < 0 somewhere. In all the cases that we have investigated (numerically), 
the eigenvalue of the tunnelling negative mode remains negative. This supports the view 
that for CdL bounces the tunnelling negative mode always exists. 

Second, we will confirm the existence of an infinite tower of additional negative modes 
when Q becomes negative. They were already conjectured to exist by LRT in their early 
work [23] , and here we provide some explicit numerical examples. These additional negative 
modes, which are always present in the region where Q < 0, remain puzzling: on the 
one hand, it does not seem possible to remove them by a canonical transformation (see 
appendix C), but on the other hand the associated instantons do not appear to be physically 
signihcantly different from those where Q > 0 everywhere. The existence of this additional 
set of modes remains somewhat mysterious and will require further clarification. 
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II. DESCRIPTION OF FALSE VACUUM DECAY IN EUCLIDEAN APPROACH 

Let’s consider the theory of a single self-interacting scalar held which is dehned by the 
following Euclidean action 

Se = J d^x[^d^ipd>^ip + V{lp)] . ( 1 ) 

Further, let’s assume that V{(p) is an asymmetric double-well potential and that it has a 
local minimum (false vacuum) at some (p = (pj, an absolute minimum (true vacuum) at 
ip = ipt and a local maximum (top) at some ip = p>top, such that ipt < (ptop < Vf- 

The energy Eq of the lowest energy state localised around the false vacuum gets a cor¬ 
rection due to quantum tunnelling effects, 


Ef — Eo — j . ( 2 ) 

It turns out that the correction is purely imaginary 7 = i| 7 |, which is a sign of metastability 
and shows that T = I 7 I actually describes the decay width of the false vacuum. This decay 
width is given by the functional integral 

7=^1 ( 3 ) 

where V_ is a normalisation factor. In the quasi-classical approximation the functional 
integral Eq. (3) can be evaluated by considering small perturbations about the classical 
saddle point Euclidean solution known as the “bounce”. The Euclidean action can be 
expanded as 

Se = + S^^'>{6ip)\^=^, . (4) 

The normalization factor is the same functional integral calculated about the false vacuum 
ip = ipf. So, for 7 we hnd the Arrhenius formula 

7 = Ae-'^ , (5) 


with 




( 6 ) 
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and A being ratio of the corresponding integrals: 

The qnadratic action for the 0(4) symmetric conhgnrations takes the form 

= 27r^ j 7]^dr]S(pdF6(p , with Op = • (§) 

So, the mode eqnation diagonalising the qnadratic action Eq. (8) has the form of a 
Schrodinger eqnation 

- -/ + = XJ^n . (9) 

dr]^ 7] drj 

Since any pertnrbation with proper bonndary conditions can be decomposed into a complete 
set of fnnctions of the flnctnation operator O^, 


5(p ^ ^ c^dipn , 

n 


( 10 ) 


integration over Sip in Eq. (7) can be replaced by integration over Cn- 
integrals, one obtains the prodnct of eigenvalnes, i.e. determinants of 
operators 


A 




(UA 



Sdet'[-d^ + V"{p>^)]\ ^ 
dvr^ \^(iet[—+ V"{ipf)] ) 


Taking Ganssian 
the corresponding 


( 11 ) 


Here a prime indicates that the zero-modes must be omitted - as described by Callan and 
Coleman [20], a proper treatment of the zero-modes results in the H^/dvr^ pre-factors. While 
hnding bounce solutions and calculating the exponential factor in the thin-wall approxima¬ 
tion (or numerically) is a relatively easy task, the calculation of the pre-exponential factor 
in held theory is considerably more involved [32, 39, 40] - in particular, one has to deal with 
possible one-loop divergences. 
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III. BOUNCE SOLUTIONS WITH GRAVITY 


Let’s consider the theory of a single scalar field minimally coupled to gravity, which is 
dehned by the following Euclidean action 

Se = J + V , (12) 

where k = SttGn is the reduced Newton’s gravitational constant. The most general 0(4) 
invariant metric is parametrised as 


= N‘^{7])d7f + p‘^{ri)dfll , 


(13) 


where N{r]) is the Lapse function, p{r]) is the scale factor and is metric of the unit 
three-sphere, 

dfll = dx^ + sin^x(d 6 '^ -|- sin^(6')d0^) . (14) 

For the metric in Eq. (13) the curvature scalar looks like 

6 6p QpN 

p2 p2^2 p^2 + piV3 ’ ^ ^ 

where ' = d/dp. Using the ansatz (13) and assuming that p = p{rj), we get the reduced 
action in the form 


Se = Se{p,N,p) = 27r^ 



+ p^NV{p) 


3 ^ 

K kN kN 


nm ) 


(16) 


In proper-time gauge, V = 1, the corresponding held equations are 


p dV 

p + 3-p = — , 

P dp 

(17) 

= -y( 0 ^ + ''^(</ 2 )) , 

(18) 

^IT^(^-U). 

3^2 ’ 

(19) 


Now let’s assume that the potential V{p) has two non-degenerate local minima a.t p = pt and 
p = pf, with V{pi) > V{pt), and a local maximum for some p = (^top, with p^ < (^top < Pi- 
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The Euclidean solution describing vacuum decay - the bounce - satisfies these equations, 
and when V ( 99 ) > 0 one has the boundary conditions 

m = ^o, 0(0) = 0, p(0) = 0, p(0) = l (20) 

at ?7 = 0 and 

^{j]max) ~ ^rri) '^{j]max) ~ 0 , p{j]max) ~ 0 , p{j]max) ~ 1 (^ 1 ) 


at some p = pmax- This assumes the following Taylor series as 77 —)■ 0: 


p{p) 


, l/'(^o) 2 , ^ , 

+ (<^o) + 


2kV{(Po)^ 4 V'{ipo) 


V + 




192 ' 3 ' 829440 

{ipof + lQ2KV'{(pQf‘ + 180fi;E((/)o)E''(99o) + 112fi;^E(</9o)^]h^ + 0{p^) , (22) 

P - ^VMp^ - ^[^^'(<^ 0 )^ - '^V{<pof]p^ 


K 


2177280 


[405E'((^o)V"((^o) - UKV{^o)V'{^of + 16«V((^o)%'' + 0{p^) ,(23) 


where V'{(po) = etc. Similar power-law behavior is valid for non-singular bounces 

for a; —)■ 0, where x = Pmax — p- These Taylor series are required when numerically solving for 
instanton solutions, as one cannot directly integrate from 77 = 0 (given the above boundary 
conditions), but rather has to start the integration at some small value 77 = e 1 . 


IV. NEGATIVE MODE PROBLEM IN HAMILTONIAN APPROACH 

Let’s expand the metric and the scalar field over a 0(4)—symmetric background as follows: 

ds^ = {1 + 2 A{p))dp^ + p{pf {I-2^!{p))dnl , (^ = <^( 77 ) + $( 77 ) , (24) 

where a and ip are the background field values and A, 4/ and <h are small perturbations. In 
what follows, we will be interested in the lowest (purely 77 -dependent, ‘homogeneous’) modes 
and consider only scalar metric perturbations. Expanding the total action to second order 
in perturbations and using the backgronnd equations of motion, we find 




( 25 ) 
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where is the action of the background solution and 'h, $] is the quadratic action 

given below. 

The quadratic action about CdL bounces was first derived in [23] using the Lagrangian 
approach and in particular it was noted that when gravity is taken into account the corre¬ 
sponding operator Oq in front of the kinetic term contains a factor (see the Appendix for 
more details) 


Kp^V{^) .2 

Qlrt — t - 5 - — p - 7 — 

3 D 


(26) 


This factor typically becomes negative somewhere along a bounce solution. Later, using a 
Hamiltonian approach in the context of the theory of constrained dynamical systems, the 
quadratic action was rederived [27] and it was shown to have a similar structure, but with 
a different factor in front of the kinetic term of the perturbations. 


Qklt = Q 


6 


(27) 


The original quadratic action [A, T, <h] is degenerate and describes a constrained dy¬ 
namical system. Applying Dirac’s formalism as in [27], one gets an unconstrained quadratic 
action for a single physical degree of freedom (which for this gauge fixing procedure is simply 
the scalar held perturbation <h) 

Sg'm =x=y'*,4(-T(dMT) + p=(,,)(/|v>(>,).M»)l)4 . (28) 

where the factor Q was given above in Eq. (27) and the potential U is expressed in terms of 
the bounce solution as 


U[piv),piv)] = ^ + P^y'^i.v) - ^PP^V'i^)) ■ (29) 

The exact form of the huctuation operator depends on the choice of a weight function, which 
can be specihed by dehning a norm. In the context of general relativity the natural choice 


is [32] 


|$|P = J d'^x y/g = 2n'^ j 


dr] p{r])^ . 


(30) 
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The fluctuation equation diagonalizing the quadratic action Eq. (28) then has the form 

1 d ,P^ d^„ X r / N / M , /X 

where and are eigenfunctions and eigenvalues of our Dirichlet boundary value problem. 
Note that this equation correctly reduces to the flat-space equation (9) in the limit k —)■ 0, 

p^rj. 

The potential U, close to ?7 = 0, behaves as 

U = Uo + U2V^ + 0[n^], (32) 

with 

Uo = l/"((po) , and U 2 = + ^E'(<^o)^'"(<do) • (33) 

D 8 

Regular solutions (eigenfunctions) close to p = 0 then behave as 

4 = A„{l + i(V'(^„)-A)y + ^l3(V'(^„)-A)" + 2K(V'(^„)-A)l/(v><,) 

+ ri‘ + Oh']} , (34) 

with Aq being a normalisation constant. Obviously, the potential U and the regular branch 
of the wavefunctions have the same behaviour in powers of x = rjmax — p Sit the end of 
the interval, close to pmax- Whereas in general there exists also a singular branch behaving 
ag ^sing ^ some function / having the property that /(O) = 0, one can adjust 

the value of A such that the singular branch is suppressed. In fact, in this way one may 
determine the (quantised) energy eigenvalues of the eigenfunctions. 


V. REGULARITY OF THE PERTURBATIONS WHEN Q PASSES THROUGH 

A ZERO 


As was discussed in Sect. IV, the equation for linear perturbations in the Hamiltonian 
approach has the form 


1 




— + U^ 

dp 


Ad), 


(35) 
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with the potential U given in Eq. (29). Note that the function Q tends to 1 at the ends 
of the interval [0,77maa;]; but for some bounces it can become negative for some interval of 
rj. The case where Q becomes negative has long been considered as puzzling [27, 29, 32] 
because Eq. (35) then shows an apparent singularity when Q passes through zero. We will 
now show that when Q{fj) = 0 the general solution of (35) is nevertheless across the point 
fj. We treat the cases where Q{fj) 7 ^ 0 and Q{fj) =0 separately. We will assume that Q and 
the other background quantities are at least functions of rj. 


A. Qii]) 7 ^ 0 case 


Let X = 7] — rj. We can re-write (35) as 


Q4 + ( Q-3^Q I t+i/i = AQ'4 , 


(36) 


with 


- , O hnpp^V .2.2 


U = Q^U = QV" + 2nQip^ + 


+ 2kp (f 


(37) 


3 3 

We start with the consideration of the zero mode equation, i.e. we assume A = 0 in Eq. (36). 
Since we are interested in showing the regularity of the wave function in the vicinity of a: = 0, 
we re-write (36) as 

(38) 


where we denote 


xA(x) + Bix) <h H—C'(x)<h = 0 , 

q 


q = Q{p), q = Q{p),... , 
u = U{f]), u=U{f]),... , 

A{x) = 1 + 0{x) , 

B{x) = l + bx + Olx"^), 6=4 — 3-, 


q oP 


q p 


u 


Cix) = l + cx + C>(x^), c= —. 

u 


(39) 

(40) 

(41) 

(42) 

(43) 


We tentatively solve the equation by keeping the dominant behaviour of the coefficients first. 


u 

a:<h + <h + -<h = 0. 

q 


(44) 
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The general solution of this equation can be written in terms of Bessel functions. The leading 
behaviour of the solution is of two different kinds: 

/ ^ \ 

<h = a 1- X -In |x| + ... + /3x^(l + ...) . (45) 

V g 2g2 J 

However, the second sub-leading term in the hrst asymptotic solution is not relevant, since 
it contributes a term of order x to the perturbation equation, which is the same as the 
dominant order we neglected in going from (36) to (44). To hnd the right expression for the 
sub-leading term and make sure that is well defined at a: = 0 , we set 


$ = a 



/ 

—X + X'lp 

Q 


The asymptotic equation for ^|J reads 


(46) 


.. , u / u \ 

— x'^'ip — x'lfj + 'il) ^— { c — b -x = 0. (47) 

Q \ qj 

Discarding the solution ~ 1/a:, which simply changes the value of a, the leading behaviour 
of the solution is the same as in (45), albeit with a different coefficient: 


$ = a ^1 — ^a: + ^ — 5 — a:^ In |a:| + .. .^ + /9a:^(l + ...) . (48) 

However, a closer analysis of the non-homogeneous piece appearing in (47) shows that its 
coefficient vanishes due to the background field equations. This means that one should move 
to the next-to leading order for the non-homogeneous piece, and the solution for 'ip is of 
order x with no logarithmic corrections: 


$ = a ^1 — ^a: + 0{x‘^)j + /da:^(l + ...). (49) 

This completes our proof: the apparent singularity in the perturbation equation forces a 
fixed relation between $( 77 ) and $( 77 ), but does not cause any divergence in the solution, at 
least for $ and its hrst two derivatives. 

We now note that, if one considers the eigenvalue equation (35) with A 7 ^ 0 the result 
above is still valid, since this merely amounts to a change in C'(a:) by a term of order 
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Therefore, the generic solution of the eigenvalue equation is also non-singular ai rj = rj. 


B. Q(t?) = 0 case 


In this case, we can’t simply borrow the previous results, as g = 0. The held equations 
alone do not determine the leading degree of the Q polynomial around rj = fj. Keeping it 
general, we will write 

0 = + 0(x*+=), k>l. (50) 

We require fc > 1 to be an integer. At hrst look, one would think that U = (Q'^U) should 
not generally vanish at rj = f], so that the leading terms of the perturbation equation would 
read 

• ?/ 

+ + = 0 , (51) 

where again u = U{fj) is generally non-zero. This would be a problem, as the equation 
above admits solutions which diverge near r] = fj. A closer analysis shows that, in fact, the 
behaviour of U near f] is related to that of Q itself. Let us show that, under the condition 
(50), one has 

• (52) 

We have 


U = Q {V + 2Kif^) + u , 


U = 


bKppifV ,2.2 

-—--|- zKjP (p 


(53) 

(54) 


The hrst two terms in (53) separately satisfy the bound (52). Concerning the last piece, 
using the background held equations one can show that 


u = 


pQ 3Q^ 


P 


Kp^lp'^ 


(55) 
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This relation is valid everywhere, not only at rj = rj. Taking its derivatives up to order k — 2, 
one easily deduces that (remembering that ip{fj) 7 ^ 0 since Qifj) = 0) 

MW(r^) = 0: n = 0,...,A:-2 , (56) 


and hence (52), 


U = 


u 


(k-i) 




X 


k-1 


+ 


U 


(k) 


k\ 




Now, dividing through by x^ ^ and keeping the leading terms near r] 
equation very similar to Eq. (44), namely 


(57) 


rj, we obtain an 


— X + Zed) + k —77^<h = 0 . 

qik) 


(58) 


The asymptotic form of the general solution can be obtained following the previous proce¬ 
dure, Ending: 

$ = a ^ - qO^x + + ■■ -j -k fix {1 + ...) , (59) 

/ f.ik-l) Jk+l) d\ 7 /(^) 

d = ^ k^^ + - 3^ - „ , ,,, . (60) 

Note that in this case (i.e. for k > 1), the In |x| term right away does not appear in the 
expansion.^ We conclude that also in this case, the apparent singularity in the perturbation 
equation does not lead to a singularity in the general solution. 


VI. NUMERICAL RESULTS 


We have studied the existence and the properties of fluctuation modes about CdL bounces 
numerically for a range of representative potentials. We have mostly used the same potentials 
as those used by Lee and Weinberg in their related study [36], as this allows us to better 
compare our results with theirs (see appendix A). 

^ We also note that the term linear in x may be absent if ^ 55 ^ shows that this can happen 

if I = 0 , which in turn can only happen when V(fj) = 0 , as can be seen by combining the expression for 
Q and the scalar field equation of motion. 
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(a) K = 0.055 


(b) «: = 0.09 


Figure 1; Profile of the bounce: (p{rj) and p{ri) for the 


potential (61) and two different values of k. 


In order to find a bounce solution we have to integrate Eqs. (17, 18) with the initial 
conditions (20). But since Eq. (17) has a regular singular point at r/ = 0 we cannot start 
the numerical integration there. Instead we start the numerical integration at some small 
f} = e, with initial conditions at this point provided by the Taylor series Eqs. (22), (23). This 
is a one-dimensional shooting problem with the shooting parameter being the initial value 
of the scalar held, (po- For the compact bounces under consideration, one has to choose ipo 
such that at ?7 = rjmax one gets the behaviour given by Eqs. (21). Since the initial conditions 
are given with hnite precision, at some rj < rfmax any such numerical solutions will start 
to deviate from the exact solution, because of a mixture with the singular branch. For 
some potentials, one needs to adjust the initial values to very high precision to construct 
the bounce solution. The presented Taylor series of Eqs. (22)-(23) guarantee precision of 
about 10“^° for e = 10“®. In the following, we specify different potentials and report our 
numerical results regarding concrete examples for various potentials and varying values of 
K. The numerical method that we employed is the Runge-Kutta of fourth order with an 
adaptive step size for the background quantities and with a hxed step size for the mode 
functions, and we implemented it in the C-I--I- programming language. 
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Figure 2: The kinetic pre-factor Q{ri) for the potential (61) and different values of k. 
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Figure 3: Lowest eigenvalues Aq of the fluctuation equation (31) versus n for the potential (61). 

A. Example 1 

As a first example, we have chosen the double-well potential also studied in [36], namely 

c {(p) = (v? - 3) ^ 

In this potential the true vacuum at <pt = 0 and the false vacuum at ~ 2.81 are separated 
by a potential barrier with local maximum at iptop ~ 1.69. 

We are interested in bounce solutions in this potential, but for different values of the 
gravitational coupling k. Fig. 1(a) shows the CdL bounce solution for this potential when 
K = 0.055, while Fig. 1(b) shows the bounce solution for k = 0.09. It was claimed in [36] 
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that a single tunnelling negative mode is present for small values of k such as k = 0.055, 
but that no negative mode is present for larger values, namely for k = 0.057, k = 0.07 or 
K, = 0.09. By contrast, we hnd that a single tunnelling negative mode exists in each of these 
cases. Note that the quantity Q is always positive for these choices of k, see Fig. 2, and 
that the (negative) eigenvalue of the lowest mode evolves continuously as a function of the 
gravitational strength k, as shown in Fig. 3. Thus, for this example, the story unfolds just 
as in the case of tunnelling in the absence of gravity. 

B. Example 2 

As the second example we choose the potential (cf. [36]) 

= + (62) 

This potential has two local minima separated by a barrier, as long as the parameter B > Bcr, 
with Bcr ~ 0.52. A CdL bounce exists for B > Bhm, where Bhm ~ 0.55. Varying the 
parameter B is instructive because for the potential (62), the factor Q is always positive along 
the bounce trajectory in the case of small values of B, while for B > Bq with Bq ^ 3.22, the 
bounce solution develops a region of negative Q. Therefore this potential is ideally suited 
for studying the case of principal interest here, namely what happens when Q becomes 
negative. We analysed the solutions for selected parameter values B varying between 1 and 
10. Our main hnding is that nothing special happens with the tunnelling negative mode 
when Q becomes negative (as can now be expected, given the analytic treatment in section 
V). Fig. 4(a) shows the bounce prohle for B = 3, while Fig. 4(b) displays the corresponding 
bounce when B = 10. Note that factor Q{ri) is always positive for B < Bq and is negative 
in some interval of r] for B > Bq, see Fig. 5. In this setup we investigated the tunneling 
negative mode and the hrst few positive modes in detail. Fig. 6 illustrates the dependence 
of the eigenvalue of the tunnelling negative mode i/jq as well as the hrst two excited positive 
modes "01^2 on the value of B. The characteristic prohles of these modes are shown in Fig. 7 
for the particular value of B = 10. Note that the eigenvalues evolve continuously as a 
function of B and in particular that the eigenvalue of the tunnelling negative mode always 
remains negative, regardless of whether Q is positive everywhere or negative somewhere. 
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(a) B = 3 


(b) B = IQ 


Figure 4: Profile of the bounce: 


ip{r]) and p{r]) for the potential (62) and different choices of B. 





{&) B = 2 and B = 3 


(b) B = 3.22 (Q < 0 briefly) 


{c) B = 5 and B = 10 


Figure 5: The kinetic pre-factor Q{rj) for the potential (62) and different values of B. For B ^ 3.22 
Q becomes negative along some interval of Euclidean time rj. 


For values of i? > So, the kinetic pre-factor Q has negative values in some interval and 
infinitely many additional negative modes are expected to have support here [23, 36]. For 
the hrst time, we will exhibit such modes explicitly here. The wave functions of the hrst 
two additional negative modes are shown in Fig. 8 for the particular value of S = 10. 
As expected, these modes have support in the region over which Q is negative. The two 
modes shown here have eigenvalues A_i = —18.22 and A _2 = —78.45 respectively. As is 
intuitively clear, the magnitude of these eigenvalues is inversely proportional to the length 
of the interval over which Q becomes negative. In Fig. 9 we show how the eigenvalues of 
the hrst two additional negative modes evolve as the parameter B in the potential is varied. 
The eigenvalues already reach A_i = —907, A _2 = —9537 at B = 3.4 and tend to —oo as 
B approaches Bo from above, while for B < Bq these additional negative modes disappear 
altogether as Q remains positive throughout. Note that over a substantial range of B the 
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(b) Ao 


Figure 6; The first two positive eigenvalues Ai and A 2 and the first (tunnelling) negative eigenvalue 
Ao of the fluctuation equation (31) versus B for the potential (62). 




Figure 7: The wave functions of the first two excited positive modes and of the tunnelling negative 
mode in the case of 5 = 10 for the potential (62). Note the separate axes for the two positive 
modes on the left and right of graph (a). 


magnitude of the eigenvalues of the additional negative modes is vastly larger than those of 
the tunnelling negative mode and the hrst couple of positive modes. 

We also studied the dependence of the above cases on the gravitational strength k. As an 
example, we hxed B = 3 and decreased the value of k. Fig. 10 shows plots of Q for different 
values of k when B = 3. Fig. 11 shows the dependence of Aq with respect to different 
choices of k. Once again, we can see that nothing special happens to the tunnelling negative 
mode as Q reaches negative values, and in particular the eigenvalue evolves smoothly. As 
before, we would expect the existence of an additional tower of negative modes wherever Q 
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Figure 8: The wave functions of the first two extra negative modes in the case of -B = 10 for the 
potential (62). 
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(a) A_i and A _2 

Figure 9: The first two extra negative eigenvalues A_i and A _2 of the fluctuation equation (31) 
versus B for the potential (62). 



is negative. 


C. Symmetric potential 

As the third example we consider a symmetric potential. This case is interesting because, 
in flat space, hnite-action Euclidean solutions in a symmetric potential have at most zero 
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(a) K = 0.9 (b) K = 0.7 (c) k = 0.5 


Figure 10: The kinetic pre-factor Q{r]) (plotted only over the region of greatest interest) for de¬ 
creasing values of n with fixed parameter B = 3 for the potential (62). 


- 0.15 


Figure 11: Dependence of the eigenvalue of tunnelling negative mode Aq on different choices of k 
for B = 3 for the potential (62). The eigenvalue evolves smoothly and remains negative. 

modes. When gravity is included, it was first found in [41] that the corresponding instantons 
develop a negative mode. For the numerical investigation, we choose the potential 

V(^) = B(,^^ - j)' + (63) 

and we vary the parameter B in the interval [2,4]. Once again we hnd a single tunnelling 
negative mode in all cases, as well as a tower of additional negative modes for parameter 
values for which Q becomes negative. Fig. 12 shows the kinetic pre-factor Q for selected 
values of B and Fig. 13 illustrates the dependence of Aq on B. The eigenvalue evolves 
continuously and remains negative as Q hrst reaches zero and then takes on negative values 
along an interval of Euclidean time r] along the instanton. 
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(a) B = 2 and B = 2.5 (b) B = 2.75 and B = 3 

Figure 12: The quantity Q{ri) for different values of B for the symmetric potential (63). 


- 0.13 


- 0.135 


- 0.14 
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Figure 13: Dependence of Aq on different choices of B for the symmetric potential (63). 



VII. DISCUSSION 


We investigated tunnelling transitions with gravity with emphasis on the cases where 
the factor Q, which sits in front of the kinetic term of fluctuations about the instantons, is 
negative in some interval of Euclidean time. We have shown that the perturbation equations 
remain well-dehned when Q crosses zero, and that the eigenvalues of existing mode functions 
vary continuously as we pass from situations with Q > 0 everywhere to Q < 0 somewhere. 
Thus, in particular, all bounces contain a single tunnelling negative mode in their spectrum 
of perturbations. We should point out that these results are partly at odds with a recent 
study of Lee and Weinberg [36, 37] - we provide a comparison of the two approaches in 
appendices A and B. 
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As we have discussed, when Q < 0 an additional inhnite tower of negative modes appears, 
in analogy with the hrst description of the problem by LRT [23]. In the present study we 
found the hrst few of these additional negative modes for concrete potentials numerically 
and determined their eigenvalues. We observed that over a substantial range of parameters 
the magnitude of the eigenvalues of the additional negative modes is vastly larger than those 
of the tunnelling negative mode. In our opinion this is another manifestation of the fact that 
they have a different nature. Furthermore, it is interesting to note that nodal theorems [42] 
are not fulhlled in this case, namely the number of nodes of the zero energy wave function 
does not count the number of negative energy states in a given potential anymore. 

As is known in certain cases, a proper choice of variables (obtained with the help of 
canonical transformations in the Hamiltonian formalism) can signihcantly simplify the de¬ 
scription of a system. In appendix C we study the influence of canonical transformations 
on these additional negative modes and provide arguments that they cannot be removed 
by a canonical transformation of variables. The consideration of canonical transformations 
shows that, even though the pre-factor of the kinetic term cannot be made positive in gen¬ 
eral, the classes of instantons for which it is negative can be changed - this suggests that 
the negativity of the pre-factor is a technical problem, not related to a fundamental physical 
cause. The existence and signihcance of the additional negative modes remain somewhat 
mysterious, and further work will be required to fully elucidate the negative mode problem. 
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Appendix A: Comparison to the work of Lee and Weinberg 

It may be useful to compare our results with those of Lee and Weinberg (LW) [36, 37], who 
recently performed a related study of the negative mode problem. Their work was performed 
in the Lagrangian formalism - in appendix B, we discuss the Lagrangian approach in more 
detail. The main difference between our results and those of LW is that there are cases 
where LW claim that no tunnelling negative mode is present. In particular, for instantons 
where the wall of the instanton is close to the maximum of p (and where the initial and hnal 
vacua are nearly degenerate), they hnd that the tunnelling negative mode is absent. We 
disagree with these results, as we have explicitly demonstrated the existence of the tunnelling 
negative mode for such cases, even for the case of exactly degenerate initial and hnal vacua (a 
possible explanation of the discrepancy between our results and those of LW is discussed in 
appendix B). Moreover, we have provided a proof that solutions to the eigenvalue equation 
evolve smoothly from cases where Q is everywhere positive to Q going negative somewhere. 
This analytic result boosts our conhdence that our results are correct. We note that LW also 
conhrm the existence of an additional inhnite tower of negative modes in regions where Q 
is negative, though it should be emphasised that in the Lagrangian approach it is the factor 
Qlrt = fhat is the relevant one, rather than the function Q = 1 — that 

appears in the Hamiltonian formalism used here. This makes a difference, as Qlrt becomes 
negative near p = 0 for all bounces (and thus all bounces will have an inhnite number of 
negative modes in the approach of LW), while the function Q is only negative for certain 
classes of bounces. We certainly do agree with LW that these additional modes are not 
properly understood yet, and that their nature and physical signihcance must be studied 
further. 


Appendix B: Negative mode problem in Lagrangian approach 


We will briehy describe the derivation of the equations for linear perturbations in the 
Lagrangian approach. For the 0(4)—symmetric situation, Eq. (24), after hxing the gauge 
T = 0 and using the constraint equation to eliminate A, one gets the quadratic action [23] 


= 27r^ 


p^drj 

‘^Qlrt 


^.2^2 _ + {QlrtV" + ^«pV'2)4.2] . (Bl) 


0 ( 2 ) 


3 
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Integrating by parts and using the background equations of motion, one obtains 

sf = j (b2) 

with the potential 

_ p^V" fi:p2p2^'2 ^pp^V' 

Qlrt ^Qlrt ‘^Q\rt 

The exact form of the fluctuation operator depends on the choice 
can be specihed by dehning the norm. With the choice 

d'^x y/g 4)^ = 27r^ j dp p{p)^ , (B4) 

the fluctuation equation diagonalising the quadratic action Eq. (B2) has the form 

-, (B5) 

p-^dp Qlrt dp 

where and A„ are eigenfunctions and eigenvalues of the Dirichlet boundary value problem. 
The potential f/$ close to p = 0 behaves as 

with 

I'f = V"(p„) - , and £;P> = + \v' (p„)V'"(p„) . (B7) 

Regular solutions (eigenfunctions) close to r/ = 0 then admit the asymptotic behaviour 

^ = -Do{l + — ^) {pQ — \Y + 2k{V''{pq) — A)R(93o) 

+ 1^kV'\po) + W{po)V"\po)] h" + 0[p^]} , (B8) 

with Dq being a normalisation constant. 

This gauge-hxed approach can easily be promoted to the gauge-invariant approach 
adopted by Lee and Weinberg (LW) [36]. For the gauge-invariant variable 



(B3) 


of a weight function which 


X = d* + pp'i 


(B9) 
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one gets a simple quadratic action 

Sf = 2^^ I f/dri + lu^X^) . (BIO) 


with the potential 


^ V” KpifV 2KpP‘ KppV 

^ Qlrt ^QIrt ^pQlrt ^Qlrt pQlrt 


4:K,V pQlrt 
^Qlrt pQlrt 


Choosing the norm to be given by 


(Bll) 


I Ixl |2 = j d^x ^ = 27r2 j dp p{pf , (B12) 

the fluctuation equation takes the form 

“ + U^Mp), p{p)]xn = KXn , (B13) 

P'^dp Qlrt dp 

where Xn and are eigenfunctions and eigenvalues of the Dirichlet boundary value problem. 
For completeness, we will also write out the asymptotic behaviour: the potential close 
to 7] = 0, behaves as 

+ 0{p^) , (B14) 


with 


up = V'(^„) - , 

K K 2/^^ 1 

up = + -V(^„)l/'"(43„) . (B15) 

Regular solutions (eigenfunctions) close to ?7 = 0 behave as 

X = Xo{l + - ^K-V{po) - \) p^ + Y^[9V'{po)V"'{po) + 9iV"{(po) - 

- 6kV{po){5V"{po) - 9 A ) + 3kV'^{po) + 8k^V\po)] P^ + 0[p^]} , ( B 16 ) 

with Xo a normalisation constant. 

We will finish with a simple, but important, remark: the definition of y is gauge invariant. 
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but let’s pick the typical gauge where one sets tl' = 0, so that we have the relation y = p<h. 
Now when describes the lowest eigenmode, i.e. the negative mode, it has a dependence 
on Tj in the shape of a typical bell curve, without nodes. However, note that at the middle 
of the instanton p = 0, so that y will thus necessarily possess an extra node compared to 
<h, even for the lowest-lying eigenfunction! The lowest-lying eigenfunction is peaked near 
the wall of the instanton, i.e. near the region where the scalar held p varies the most. If 
this region is well-separated from p = 0, one may therefore not see the node, and it may 
look like x is nodeless. However, when the wall is located near p = 0, the negative mode, 
expressed in terms of x, will have one node and will naively look like the hrst excited mode. 
In this case, one must be careful not to discard it erroneously. This observation may in part 
explain the discrepancy between our results and the claims in [36]. 


Appendix C: Canonical Transformations 

The classical equations of motion of a system can be derived by requiring the stationarity 
of the action, the integral of the Lagrangian function. 


55' = 0 —)■ (5 / {pq^t — H) dt, 


(Cl) 


where {q,p) are the canonical coordinate and momentum, while H is the Hamiltonian of the 
system. The resulting equations of motion are 


dH 


dH 


o 5 P,t O 

op oq 


(C2) 


We could equally well describe this system with different canonically conjugate coordinates 
{q,p) and a new Hamiltonian K, 


e f r~~ ~ dK ^ dK 

0 J {pq,t - K)dt, q^t = = 


(C3) 


These two descriptions are equivalent if the two variational principles agree, i.e. if 


dF 

pq,t -H = pq^t -K + —, 


(C4) 
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where it is important to note that one may add an arbitrary total time derivative to the 
right hand side of the above equation (we are ignoring scaling transformations here, which 
would have arisen by allowing the right hand side to also be multiplied by an overall factor). 
This function F may depend on various combinations of the old and new coordinates - for 
us, the particular choice 

F = q-p + F 3 {p,q,t), (C5) 

will be useful, where F^ denotes an arbitrary function of the old momentum, the new coor¬ 
dinate and of time. With this choice for F, Eq. (C4) becomes 

dF^ dF-, dF-, 

- H = pq^t - K + p^tq + + -^P,t + ■ (C6) 

Since we are treating the old and new coordinates and momenta as independent variables, 
we can see that the equation above is satisfied if we identify 


q = 


dp ’ 


p = 


dq 


K = H + 


dF. 


dt 


(C7) 


Now let us specialise further to the choice 


F3 





(C8) 


where c{t),d{t) are explicit, a priori arbitrary, functions of time. We have included a time- 
derivative term in the denominator of the second term such that c, d are dimensionless 
functions of time. This leads to the canonical transformations 


<1= q= dq - Fp 

p = dp P = dP 


K = H{q,p) 


ld{d/cd) ^2 
2 dt ^ 


(C9) 


The important change is in the last term: the momentum-squared part of the Hamiltonian 
(which is responsible for the kinetic term in the Lagrangian description) gets modihed by 
this transformation. 

Now let us look in more detail at the case we are interested in. The Euclidean action Se 
was given in Eq. (28) (where we will drop the overall factor of 27r^ arising from the volume 










of S^). Its Lorentzian equivalent is given by 



(CIO) 


where Q, U are here understood to be given in terms of ordinary time t = —irj, e.g. Q = 1 + 
p^(f)\/Q. With the canonical momentum given by 11 = the corresponding Hamiltonian 

reads 

H = -Ct = . (cii) 

If we now apply the canonical transformation above, with (<h, H) = (g,p) —)■ {q,p), we obtain 
the new Hamiltonian 


K = - 

2 


d‘^Q 




H—— [d/c^t),t 


P H- jpq + 


c,td 


2 d2 


(C12) 


This theory can equally well be described by a Lagrangian of the form 


C = 






(C13) 


where the explicit form of the mass term is not important for our present purposes. Trans¬ 
forming to Euclidean time leads to the Lagrangian 


Ce = 


{d/c),. 


$2 _ 


$2 


(C14) 


where now Q,U take their Euclidean expressions (27), (29). Note the most important 
consequence of the canonical transformation: the coefficient in front of the kinetic term has 
changed. Can it be made positive? Unfortunately, this seems impossible in general: in order 
to make the kinetic coefficient non-singular, we should take 1/c oc Q (with no other zeros 
than those of Q) in order to cancel the divergence in U - see Eq. (37). For some cases of 
interest, e.g. the potential in Eq. (62) with 5 = 4, Q^U is positive everywhere. Then, in 
the interval where Q < 0, the second term in the kinetic coefficient will make the coefficient 
even more negative. Hence the sign of the overall coefficient will depend on the last term, 
— {d/c)^ri- At the moments where Q = 0, this last term will contribute in proportion to —dQ. 
Keeping in mind that Q has different signs at both ends of the Q < 0 interval, it implies that 
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d must change sign, and hence pass through zero during that interval. But when d passes 
through zero the pq and terms in Eq. (C12) will be divergent, rendering the new theory 
singular. There is nevertheless one lesson that one can draw from this discussion: even 
though this type of canonical transformation cannot make the kinetic coefficient positive 
and non-singular, it can change the class of instantons for which it is negative. For instance, 
for instantons where Q > 0 throughout, one may transform to a theory where the coefficient 
is negative somewhere. This suggests that the sign of Q might be “just” a technical issue, 
and not something with fundamental physical signihcance. 

There is another set of canonical transformations worth considering. These arise by 
considering the generating function 


F = Fi{q,q,t ), 


(CIS) 


leading to 


P = 


dFi 

dq 


P = 


dFi 
dq ’ 


K = H + 


dFi 


If we choose Fi = f{t)qq, we obtain 


(C16) 


q = 



P = f{.t)q, 


K = H+{-qP)^t = H. 


(C17) 


This transformation interchanges coordinates and momenta. Above, we saw that there 
are examples, such as the potential (62) with 5 = 4, where the potential is everywhere 
positive. Hence one may ask whether it is possible to use the transformation just described 
to exchange a theory with a (partly) negative kinetic term and a positive potential for one 
with a positive kinetic term and a (partly) negative potential (see also [29]). This would 
remove the unwanted inhnite tower of negative modes. Specialising to the model of interest, 
we hnd the resulting Euclidean theory 


Cr = 


P 


2p^U 


•9 QP 9 
2p3 


(CIS) 


Note that now the kinetic term is indeed manifestly positive when 17 > 0 ! However, as 
it stands the theory is still singular, due to the singularities in U. We can remove those by 
choosing / oc l/Q, with no other poles than those of 1/Q. This renders the kinetic term 
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regular, but unfortunately leads to a new potential oc l/Q, which is divergent when Q passes 
through zero. Thus, even though we can cure the kinetic term by this procedure, we are 
in effect shifting the singularity to the potential. Thus, unfortunately it appears that it is 
impossible to cure the negative Q complications by a canonical transformation. 
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